Tunnelling of condensate magnetization in a double-well potential 
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We study quantum dynamical properties of a spin-1 atomic Bose- Einstein condensate in a double- 
well potential. Adopting a mean field theory and single spatial mode approximation, we characterize 
our model system as two coupled spins. For certain initial states, we find full magnetization oscil- 
lations between wells not accompanied by mass (or atom numbers) exchange. We identify dynamic 
regimes of collective spin variables arising from nonlinear self-interactions that are different from 
the usual Josephson oscillations. We also discuss magnetization beats and incomplete oscillations of 
collective spin variables other than the magnetization. Our study points to an alternative approach 
to observe coherent tunnelling of a condensate through a (spatial) potential barrier. 
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In different branches of physics, various systems larger 
than atomic size have been examined to have signatures 
of purely quantum effects particularly in quantum coher- 
ence and tunnelling Q, |^ y{ . The macroscopic quantum 
tunnelling of magnetization (MQTM) 0, |^ Q is such 
an effect that has been vigorously sought after and has 
been claimed to be observed in molecular nanomagnets of 
Mni2 — acetate and Fes compounds 0, ll| ■ The surpris- 
ing event of MQTM was also believed to occur in iron 
storage protein ferritin 0. Apart from its fundamen- 
tal significance at the interface of quantum and classical 
realms, MQTM in a large spin system is expected to play 
a key role in realizing technological advances for future's 
quantum computers and for the newly emerging field of 
spintronics. Being macroscopic quantum objects them- 
selves, optically trapped atomic Bose-Einstein conden- 
sates (BECs) [iO, [i3 are also proposed to be promising 
candidates for MQTM ^^1) though strictly speaking, one 
may classify such systems as mesoscopic as well. 

In optical traps hyperfine spin degrees of freedom be- 
come accessible and for a condensate of atoms in the F — 
1 hyperfine ground state, its order parameter is a three 
component {Mp = — 1,0,+1) spinor. The tunnelling as- 
sociated with the internal (spin) degree of freedom can be 
induced by a transverse magnetic field when it exceeds 
a critical value depending on the strength of magnetic 
spin dipole-dipole interaction . Such Josephson oscil- 
lations with internal degrees of freedom can be compared 
to the well-known superfluid "^He Josephson oscillations 
observed in recent years ^^l- The tunnelling associated 
with the external (spatial) center of mass degree of free- 
dom has been demonstrated with cold atoms [T^ ITslIlp 
as well as atomic condensates in an optical lattice [l7| . 
The underlying strong correlations in spin and motional 
degrees of freedom in such systems may eventually lead 
to novel applications in spintronics. Beyond the semi- 
classical macroscopic Josephson oscillation physics, two 
tunnelling coupled condensates can also display many 
particle quantum correlations when the strength of tun- 
nelling coupling is tuned relative to the typical mean field 



interaction strength inside each well [Ty ■ While this phe- 
nomenon is well understood for scalar condensates, it is 
so far unexplored for spin-1 condensates [l9l l20l |2]| . In 
this article, we propose a different and more direct mech- 
anism for observing coherent condensate Josephson dy- 
namics based on magnetization oscillations due to an ex- 
ternal (spatial) Josephson coupling of a spin-1 condensate 
in a double well potential, building upon previous stud- 
ies of a spin 1/2 condensate in a double well potential 

Our model system consists of a spin-1 condensate in 
a double-well potential. When tunnelling through the 
potential barrier is weak, the wave function effectively 
behaves as a superposition of those localized in the left 
(L) and right (R) wells. In the mean field approxima- 
tion, we write the three component order parameter as 
ipiif,t) = (j)Li{'r)^i{i) + <i)m{f)T]i{t) with i = -,0, -I- re- 
spectively for each Zeeman state \Mf). (p^i is the ground 
state spatial wave function of the z-th spin component in 
the j^-th well {iy—L,K). Validity of the mean field approx- 
imation is discussed in the conclusion part of the paper. 
As has been studied in great detail [l^H^, both spin-1 
condensates of ^^Rb and ^'^Na atoms are dominated by 
spin symmetric interaction. This leads to the wide use of 
single spatial mode approximation (SMA) 26], whereby 
the mode function is taken to be the same of all three 
components, and is itself determined only from the sym- 
metric interactions. We will adopt this approximation 
as recent studies have fully delineated its validity regime 
[IE nil • Denote atomic density as (f) , hence we 
take = ^Jn,^{r), which leads to 



^{f, t) = yM^m + VuR^m, 



(1) 



where f,^ (/^ = ^ and = ff) is the spin-1 spinor asso- 
ciated with the condensate in the ly-th well, and should 
be interpreted as a column vector, while fl a row vector 
in our notation. 

We follow the convention of Ref. 0| to split the total 
Hamiltonian for an interacting spin-1 condensate into two 
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parts H = Hs + Ha- After integrating over the spatial 
variables and adopting the SMA introduced above we get 

y=LM. ^ 

+j(f77+^e), (2) 

Ha = \Y. ^^-^^ E flF^^flFjf-^^ (3) 



v—L,R j—x,y,z 



iS/A) 



where Fj are spin-1 matrices. We have defined 
cs/A J dr'nii'r), the single well ground state energy 

= j dr (J^ [VV^'' + V^V^^ , (4) 

and a positive tunnelling coefficient 

J = j dr(^^V^-\I^, + ^V^/^ . (5) 

The spin symmetric and asymmetric coefficients are cs = 
47r/i^(ao + 2a2)/3Af and ca = 47rfi,^(a2 — ao)/3A/ in terms 
of the scattering length oq (02) in the total spin channel 
of (2) of the two colliding spin-1 atoms. The tunnelling 
coupling is the same for all spin components as the po- 
tential is assumed spin independent, i.e., Vij = VSij for 
simplicity. The terms involving small spatial overlaps be- 
tween the wave functions in each well are neglected |l9j . 

The time dependent Schrodinger equations for the 
spmors are obtained from ihdi/dt = SH{i, fJ, i^,rf)/6^, 
which gives in a compact form, 



ih 



d4 

dt 
df] 



[A). 



(6) 



the equivalent coupled Gross-Pitaevskii equation for the 
two spinors. We have defined a short hand notation h^, = 

f'l ^ f'l with /'^ = (/-, -/o, /+). Both e and 77 are 
normalized to unity, so we attempt to find the stationary 
states satisfying 



dt\v J \V 



(7) 



In this short article, we will focus on the tunnelling 
dynamics in a symmetric double well with = £r — 
e,X^f^ — — Xs,X^j^^ = = Xa- In this case, 
inter-well tunnelling couples the same spin components, 
while intra-well interactions cause spin mixing, coupling 
components |-|-) and |— ) to |0). We have performed ex- 
tensive simulations and discovered a variety of interesting 
solutions: self-trapping for each spinor component, for 
combinations of different spinor components, as well as 
general nonlinear Josephson type complete oscillations. 
To illustrate these results, we will selectively display some 




FIG. 1: The Josephson phase diagram for J — 0.02 and Xa = 
—0.01. The trajectories of Eq. @ stay on closed paths. Xa < 
(> 0) refers to a ferromagnetic (antiferromagnetic) spin-1 
condensate as for ^'^Rb (^^Na) atoms and |Aa| |As|. 



oscillation dynamics that might be useful to guide ongo- 
ing experimental efforts. For these solutions to be phys- 
ically meaningful, we have also studied their stabilities 
by a standard linearization analysis, and only stable so- 
lutions are considered. 

We now try to formulate the evolution of the most 
easily measurable quantity: magnetization in each of the 
symmetric double well. Since the total magnetization 
of our system is a conserved quantity, we therefore only 
need to consider the magnetization in the one of the wells, 
e.g. the left well. We further consider a special initial 
state r]+ — rjo — S^q, and ?7_ = The density 

matrix of the system then becomes p = ^* ^ £^ with 
Pij = for i,j = -|-,0, — . It may be noted that p 

stands for single-particle density matrix, as the many- 
body density matrix is beyond the scope of a mean field 
analysis. The dynamical variables contributing to the 
evolution of magnetization M are given by 
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[(1 - no)^ -f M)][(l - no)^ - M]/4, which comes di- 
rectly from the definitions of i?Oj lo: M, and the 
normalization condition for /. Quite generally, there 
is a conserved quantity, in this system. Further 

restricting the initial state to ^0 = we find a simple 
situation analogous to the standard Josephson junction 



for a scalar condensate in a double well |19|, |2 

M = 4J/o, 
Rq = -2XaIoM, 
io = 2XARaM - JM. 



(9) 
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FIG. 2: The phase portrait for J = 0.0051 and Xa = -0.01. 



FIG. 3: Dependence of the period to the 2J/\Xa\ 



This set of equations can be interpreted as the familiar 
two state optical Bloch equations, albeit nonlinear ones, 
with i?o and /q being the analogous real and imaginary 
parts of the atomic dipole moment. M then acts as the 
population inversion. The fixed point /o = M = is a 
center, which is not an attractor. When ,/ < \\a\, there 
exists another fixed point /q = 2A4i?o — J = 0, which is 
unstable since one of its three characteristic frequencies 
is positive. Since M is physically bounded between 1 
and —1, we anticipate there exists a limit cycle in this 
system. For ^+ = 1, we find i?o = (1 - M2)Aa/4J. The 
phase diagrams for the case of full Josephson oscillations 
is shown in Fig. ^ Alternatively, we can construct two 
dimensional phase portraits by introducing an auxiliary 

variable 9^ = tan~^ {Iq/Ro) as illustrated in Fig. ^ 

Equation © can be analytically solved in terms of the 
Jacobi elliptic functions |2^J with the oscillation period 
given by 



2i^(^)/|A^|, 



4X(%I)/(2J), 



2J<\Xa\, 
2J> |Aa|, 



(10) 



where K{.) is the complete elliptic integral of the first 
kind. The behavior of the oscillation period with re- 
spect to 2J/\Xa\ is shown in Fig. |3| As we can see, 
for small values of 2J/|A^| the gradually increasing pe- 
riod resembles that of a nonlinear pendulum oscillation, 
which exhibits self-trapping. The sel-trapping condition 
is therefore 2 J < \Xa\- 2J/\\a\ = 1 is the critical value 
corresponding to the homoclinic orbit of the equivalent 
pendulum being completely in the top position. Beyond 
that the period decreases with increasing 2J/\Xa\, as the 
equivalent pendulum assumes a librator rotation, i.e. a 
full oscillation of the magnetization. Figure 0] clearly 
illustrates these different behaviors for the initial state 
Cl0) = (1,0,0) and ^(0) = (0,0,1). 

When the population in the Zeeman state \Mf — 0) 
becomes non-zero, even very small, for instance with 
^+ = r?_ = 0.9962, = % = 0.0872, and = ??+ = 0, 
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FIG. 4: Initial state is ^(0) = (1,0,0) and ff{0) = (0,0,1). 
e = 1.0, As = 1.0, and Xa = —0.01. Left figure shows self 
trapping when 2J < |Aa|, solid line for J — 0.001, dashed line 
for J — 0.0049. Right figure shows magnetization oscillation, 
solid line for J = 0.0051, dashed line for J = 0.01. 



interaction caused self-trapping disappears, while the full 
magnetization oscillation reappears in the previously self- 
trapped regime. Yet, surprisingly, the presence of the 
nonlinearity due to atomic self-interaction remains visible 
through dynamic phenomena analogous to self-trapping 
in the collective variables of other combinations of the 
single particle density matrix elements. In particular if 
we examine the collective variable p++ — poOi we find 
that it does not always exhibit full oscillations. As an 
example, we demonstrate this effect, unique to a spinor 
condensate, in Fig. El for — poo- With the given ini- 
tial conditions full (complete) oscillations between 1 and 
occur when the spin dependent atomic self-interaction 
is turned off, as shown by the dashed-line. Inclusion of 
the self-interaction, on the other hand, inhibits the full 
oscillation, leading to a similar self-trapping effect as in 
the non-linear Josephson dynamics of a scalar condensate 
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FIG. 5: Oscillation of the collective variable p++ — poo for 
an initial state f(0) = (0.9962,0.0872,0). Other parameters 
used are e = 1.0, Ag = 1.0, Aa = —0.01 (solid line), and 
Aa = (dashed line for the case of no spin dependent atomic 
interaction) , and J = 0.001. 

in a double well. As a general rule, we find oscillations 
now become incomplete over any regular interval. The 
values of p-1-4. — poo can even become negative as a result 
of including the self-interaction. We can summarize the 
above observations based on the phase space nonlinear 
dynamics. When is non-zero and in the self trapping 
regime when 2 J < \Xa\, our model system possesses one 
attractor, which attracts trajectories of different collec- 
tive variables depending on the initial conditions. 

In the opposite regime to self-trapping when 2 J > |Aa|, 
the system possesses two distinct attractors, leading to 
a complete oscillation of magnetization between —1 and 
1. We can adjust the system parameters to shrink the 
distance between the two attractors in the phase space, 
and generate beatings in the magnetization oscillation 
with appropriately chosen initial conditions. The math- 
ematical reason is simple. Due to the presence of the 
self- interaction in the limit of 2 J > \Xa\, the orbits of 
magnetization in the phase space now become open and 
continuously circulate around two closely spaced focuses 
(of the two attractors), leading to beating patterns. 

Finally, we provide some supporting arguments for the 
mean field treatment of the tunnelling dynamics. The va- 
lidity of the mean field theory can be intuitively expected 
due to the large numbers of atom in our model. At iV ~ 
10^ atoms, our system is beyond the microscopic regime, 
and becomes at least mesoscopic if not completely macro- 
scopic. We note that the period of magnetization oscilla- 
tion, given in Eq. (^UJ, has a simple dependence on the 
number of atoms r ~ 1/A''. Following the same analysis 
as in Ref . j we can estimate the effect of atom number 
fluctuations by \[t{N ± VN) - t{N)]/t(N)\ ~ 1/VN, 
which is less than 0.003% for the present system. We 



have simply taken the standard deviation a{N) of the 
atom number fluctuation to be ^ consistent with 
the mean field theory. We can also estimate the absolute 
phase diffusion times of our system due to atom number 
fiuctuations. The spin symmetric interaction in a spin- 
1 condensate, gives rise to a diffusion time of Tc ^ 20 
(seconds) for the overall condensate phase; while the spin 
asymmetric interaction, gives rise to a diffusion time of 
~ 5000 (seconds) for the relative phases between dif- 
ferent condensate components. In the above estimates, 
we have simply taken ri'^' ~ N/[a{N)cs{n)] as for a 
scalar condensate and ri^'' w N/[a{N)cA{n.)] for the 
spin mixing dynamics p^ . We have also assumed an 
average condensate density of {n) = 1.7 x 10^'^(cm)~'^ 
for N = 10*" '^^Rb atoms, which corresponds to a typi- 
cal Josephson oscillation period of ~ 1 (second). Thus 
we conclude that macroscopic magnetization tunnelling 
in a spinor condensate as modelled here can be exam- 
ined faithfully with the mean field theory. We further 
note that in a recent experiment with a smaller number 
of atoms N ^ 1000, the nonlinear Josephson oscillations 
observed were in complete agreement with the predic- 
tions of the mean field theory 32] . 

In conclusion, we have studied the double well system 
of a spin-1 atomic condensate using the mean field the- 
ory and SMA. We have characterized the various regimes 
of the nonlinear dynamics for the resulting macroscopic 
quantum tunnelling phenomena. In addition to features 
commonly attached to the scalar Josephson dynamics of a 
single component condensate in a double well, e.g. coher- 
ent atomic population oscillation and macroscopic self- 
trapping [1^ I23 , we have found interesting effects solely 
due to the spinor nature of the condensate such as the 
macroscopic oscillation and self-trapping of condensate 
magnetization without net changes of total atom num- 
bers within each well. The physics of these correlated 
tunnelling dynamics is essentially the same as what was 
found before in a double well system of a two compo- 
nent (or spin 1/2) condensate [23, but can become sig- 
nificantly richer due to the additional freedom in atomic 
internal states. For instance, amplitude modulations in 
the magnetization oscillation become possible in a spin-1 
condensate. We have illustrated the tunnelling induced 
macroscopic magnetization oscillations with several ex- 
amples for a spin-1 condensate of ferromagnetic interac- 
tions. Our results highlight the coherent Josephson dy- 
namics of a spinor condensate in a double well, and lay 
the ground work for future studies of the quantum state 
of atoms in such a system. 
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